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Abstract 

The Liouville and first Bogoliubov hierarchy equations with derivatives of noninte- 
ger order are derived. The fractional Liouville equation is obtained from the conser- 
vation of probability to find a system in a fractional volume element. This equation 
is used to obtain Bogoliubov hierarchy and fractional kinetic equations with fractional 
derivatives. Statistical mechanics of fractional generalization of the Hamiltonian sys- 
tems is discussed. Liouville and Bogoliubov equations with fractional coordinate and 
momenta derivatives are considered as a basis to derive fractional kinetic equations. 
The Fokker-Planck-Zaslavsky equation that has fractional phase-space derivatives is 
obtained from fractional Bogoliubov equation. The linear fractional kinetic equation 
for distribution of the charged particles is considered. 



The theory of integrals and derivatives of noninteger order goes back to Leib- 
niz, Liouville, Riemann, Grunwald, and Letnikov. Fractional calculus has found 
many applications in recent studies in mechanics and physics. Fractional equa- 
tions, which have derivatives of noninteger order, are very successful in de- 
scribing anomalous kinetics, transport, and chaos. Fractional kinetics equations 
usually appear from some phenomenological models. In this paper, we suggest 
fractional equations of statistical mechanics. To obtain these equations, the 
conservation of probability to find a system in a fractional differential volume 
element of the phase-space is used. This element can be considered as a small 
part of the phase-space set with noninteger dimension. The suggested fractional 
equations of statistical mechanics are used to derive the fractional kinetics equa- 
tions. 
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1 Introduction 



Fractional equations [H [2] contain derivatives of noninteger order [3l H] . Integrals and deriva- 
tives of fractional order have found many applications in recent studies in mechanics and 
physics. In a short period of time the list of such applications becomes long. For example, 
it includes chaotic dynamics [5l[6], mechanics of fractal media [71[8l[9], quantum mechanics 
[Tot [TT] . physical kinetics [5], [121 [IS [IH [IS], plasmas physics |T6l [17], long-range dissipation 
[TSl [T9l [20] , mechanics of non-Hamiltonian systems [211 122] , theory of long-range interaction 
[231 [Ml [25] , anomalous diffusion, and transport theory [3, [261 EZ] • 

Equations, which involve derivatives or integrals of noninteger order are very successful in 
describing anomalous kinetics [5l [6l [121 [131 [H] . Usually the fractional equations in dynamics 
or kinetics appear as some phenomenological models. In [2T] the attempt to derive the 
basic statistical mechanics equations with derivatives of noninteger order have been realized. 
Unfortunately, the fractional derivatives appear only by Fourier transform of these equations 
as it realized for the Fokker-Planck-Zaslavsky equation in [9]. 

In this paper, we derive the Liouville equation with fractional derivatives with respect to 
coordinates and momenta. To obtain the fractional Liouville equation (FLE), we consider the 
conservation of probability in the fractional differential volume element. This element can be 
considered as a small part of the phase-space set with noninteger dimension. Using the FLE, 
we get a fractional generalization of the Bogoliubov hierarchy equations. These equations 
can be used to derive fractional kinetic equations [5l [6l [121 [13] • The Vlasov equation with 
derivatives of noninteger order is obtained. The Fokker-Planck-Zaslavsky equation, which 
has fractional phase-space derivatives, is derived from the FLE. The linear fractional kinetic 
equation for distribution of the charged particles is suggested. 

In Sec. 2, we obtain the Liouville equation with fractional derivatives from the conser- 
vation of probability in the fractional volume element of phase space. In Sec. 3, the first 
Bogoliubov hierarchy equation with fractional derivatives in phase space is derived. In Sec. 
4, we consider the Vlasov equation with fractional derivatives in phase space. In Sec. 5, 
the Fokker-Planck-Zaslavsky equation, which has fractional derivatives with respect to co- 
ordinates and momenta is considered. In Sec. 6, the linear fractional kinetic equation for 
distribution of the charged particles is obtained. Finally, a short conclusion is given in Sec. 
7. 
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2 Liouville equation with fractional derivatives 



A basic principle of statistical mechanics is the conservation of probability in the phase- 
space. The Liouville equation is an expression of the principle in a convenient form for the 
analysis. In this section, we derive the Liouville equation with fractional derivatives from 
the conservation of probability in a fractional volume element. 

In the phase space R'^"^ with coordinates (x^, ...,x^^) = (gi, ...,p„), we consider 

a fractional differential volume element 

d'^V = (Txi... (rx2n. (1) 
Here, dP" is a fractional differential [28] that is defined by 

2n 

dV{x) = Y,D':,J{x){dxur, (2) 
fe=i 

where D'^^ is a fractional derivative [5] of order a with respect to Xk- 

The fractional derivative has different definitions [U [3] , and exploiting any of them de- 
pends on the initial (boundary) conditions, and the specifics of the considered physical 
processes. The classical definition is the so-called Riemann-Liouville derivative [3]. Due to 
reasons, concerning the initial and boundary conditions, it is more convenient to use the 
Caputo fractional derivatives [29l[l]. Its main advantage is that the initial conditions take 
the same form as for integer-order differential equations. The Caputo derivative is defined 
by 



1 /- /( 



n). 



Z 



D'^fix) = ^D^fix) = — ^ / , ' dz, (3) 

where n — 1 < a < n, and f^"'\z) = d^f{z)/dz'^. Note that -D°^l = 0, and D^^xi = 0, where 
{k ^ I). Using ([2]), we obtain 

d'^Xk = D'^^Xkidx^r. (4) 



Then 



From 



T{(3+l-a 
where P > a > 0. Equations (jl]) and ([6]) give 



{dxkr = {Dlxk) 'd'^Xk. (5) 



Dl4 = J}^t'\y r, (6) 



{dxkr = r(2 - a) xt^d'^Xk. (7) 
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The conservation of probability for the usual phase-space volume element is expressed as 

-dV^^^ = d[p{t,x) {u,dS)]. (8) 
For the fractional volume element ([T]), 

= t;"[p(t^a;) (u,c/"S)]. (9) 

Here, p = p{t, x) is the density of probability to find the dynamical system in d"V, u = u{t, x) 
is the velocity vector field in i?^", d'^S is a surface element, and the brackets ( , ) is a scalar 
product of vectors 

2n 2n 2n 

u = ^M,efc, rf"S = ^rf°^,efc, (u,d"S) = (10) 

k=l k=l k=l 

where are the basic vectors of Cartesian coordinate system, and 

d"Sk = d^xi... d"xk-id"xk+i... d°'x2n- (11) 

The functions = Ukit,x) define Xk components of u{t,x). In the usual case (a = 1), the 
outflow of the probability in the Xk direction is 

d[puk]dSk = D^^[puk]dxkdSk = D^^[pUk]dV. (12) 

For a 7^ 1, 

d'^lpukjd^Sk = D:^[puk]{dxrd"Sk. 
Using (mD, ([1]) and (P, we get 

d"[pUk]d"Sk = D^lpuk] {Dlxky'd'^Xkd'^Sk = {D^Xky' D^lpukjd^V. (13) 

Substitution of (fT3l) into gives 

- d-V-t = d-Vj2{D:,x,)-'D:^[pu,]. (14) 

k=l 

As a result, we obtain 

rj 2n 
k=l 

where 

= {Dlx,)-'Dl=r{2-a)xt~'Dl. (16) 
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This is the Liouville equation with the derivatives of fractional order a. Equation (iTSil 
describes the probabihty conservation for the fractional volume element ([1]) of the phase 
space. 

For the coordinates (gi, ■■■,Pn), Eq- ([T5|) is 

n " re 

k=l k=l 

where \4 = Uk, and Fk = Uk+n {k = l,...,n). The functions Vk = Vk(t,q,p) are the 
components of velocity field, and = Fk{t, q,p) are the components of the force field. 
In general, 

D^lpFk] + pD'^^Fu + FkD^^p. (18) 
If Fk does not depend on p^, and Vk does not depend on g^, then Eq. (ITTl) gives 

i + E(^^D^./^ + ^^D»=0. (19) 

k=l 

For the fractional generalization of Hamiltonian system [22] , Vk and Fjt can be presented as 

V^ = D;^i7(g,p), F, = -D^^.i7(g,p), (20) 



where H{q,p) is a fractional generalization of Hamiltonian. Substitution of (l20i) into (JT9i) 
leads to 

i + E (D.".™?./' - ^IHB^P) = 0- (21) 



dt 

k=l 



We can define 



re 



{A, 5}„ = J] {UlA B^B - BIB B;A) 



k=l 

= E i^l^^ D^P^) ' (^." ^ ^P.B - D^A) . (22) 

k=l 

For a = 1, Eq. fl22p gives the Poisson brackets. Note that 

= -{5,A}„, {i,A}„ = o. 

Using (1221) . we get fl?T|) in the form 

1^ + {p, i/}. = 0. (23) 

As the result, we have the Liouville equation for fractional generalization of Hamiltonian 
systems [22]. For a = 1, Eq. fl23|) is the usual Liouville equation. 
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3 Bogoliubov equation with fractional derivatives 



Let us consider a classical system with fixed number of identical particles. Suppose that 
A;th particle is described by the generalized coordinates qks and generalized momenta pks, 
where s = l,...,m. We use the notations = (qki, ■■■,qkm) and = {pki, ■■■,Pkm)- The 
state of this system can be described by the distribution function pat in the 2m A^- dimensional 
phase space: 

P7v(q, P, t) = p(qi, pi, q^v, Ptv, t). 
The normalization condition is 

/[l,..,iV]p^(q,p,t) = 1, 

where ...,N] is the integration with respect to qi, pi, qjy, Pat. 
The fractional Liouville equation is 

^ = -E(d^.(V,p^) + D^jF,p^)), (24) 

k=l 

where is a velocity of kth particle, is the force that acts on kth particle, and 

m 

s=l 
m 

Dp.F. = {D-.^Pk)-' D-^F, = J2 {DLp^s)-' D^F,.. (26) 

s=l 

The one-particle reduced distribution function pi can be defined by 

Pi(q,P,^) = p(qi,Pi,^) = ^[2,--,A^]p7v(q,P,t), (27) 

where /[2, A^] is an integration with respect to q2, qA?, P2, Pa^- Obviously, that the 
function fl27|) satisfies the normalization condition /[l]pi(q, p, t) = 1. 

The Bogoliubov hierarchy equations [301 EH ESI ES] describe the evolution of the reduced 
distribution functions, and can be derived from the Liouville equation. To obtain the first 
Bogoliubov equation with fractional derivatives from Eq. ^IM we consider the differentiation 
of flTTj) with respect to time 

|^4fp,....^W.fp,....^,^|^. ,28, 
Using (EID and (EH]), we get 

= -/[2, ...,iV] 5^(d«JV,p;v) + ^U-PkPN)). (29) 

k=l 
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Let us consider the integration /[q/t] over for A;th particle term of Eq. (l29l) . 
/[qfc]D^jVfep^) = /[qfc] {D^ciky' D^^iV.p^) = 1(2 - a)i[q,K-'D^^iV,p^) = 

= T{a)T{2 - a)r[cik]D^^{VkPN) = T{a)T{2 - a) (v^p^v)^" = 0, (30) 

V / — oo 

where /"[q^] is a fractional integration with respect to variables q^. In Eq. fl30p . we use that 
the distribution Ptv in the limit q^ ±oo is equal to zero. It follows from the normalization 
condition. If the limit is not equal to zero, then the integration over phase space is equal to 
infinity. Similarly, we obtain 

+ 00 



/[pfe]D°^(Ffep^) ~ (FfeP^)_[ = 0. 



Then all terms in Eq. (1291) with k = 2, N are equal to zero. Therefore, Eq. (1291) are 

^ = -/[2,...,iV](D^^(Vip^.) + D^^(Fip^.)). (31) 
The first term in Eq. flHTl) can be written as 

/[2, iV]D^jVipA.) = Vi/[2, N]p^ = D^^ (V^pi). 
For the binary interactions, 

N 

Fi = F^ + 5^Fifc, (32) 

k=2 

where F^ = F'^(qi,pi,t) is the external force, and Fi^. = F(qi, pi, q^, p^, t) are the internal 
forces. Using (1521) . the second term in (IHT]) is 



N 



/[2, iV]D^^ (Fip^) = /[2, iV] D^^ (F^p^) + (F.^Pn) 



k=2 

N 



= D^^ (F^pi) + J2 Dp,/[2, N] (Fi,p^) . (33) 

k=2 

We assume that the distribution function is invariant under the permutations of identical 
particles. Then p^ is a symmetric function, and all (A^ — 1) terms in Eq. (|33|) are identical: 

N 

J]/[2,...,iV] B^^^(f,,p^)={N-1)I[2,...,N] D^^(Fi2P^). (34) 

k=2 

Using i[2, ...,N] = /[2]/[3, A^], we have 

/[2, iV] D^^ (Fi2P^) = i[2] D^^ (Fi2/[3, iV]p^) = D°^/[2]Fi2P2, (35) 



where 

P2 = p(qi, Pi, q2, P2, t) = /[3, A^]pAr(q, p, t) (36) 

is a two-particle distribution function. 
Finally, we obtain 

^ + D;(ViPi) + D;,(F^pO =/(P2), (37) 

where 

I{p2) = ~{N-l)B^J[2]FuP2 (38) 

describes a velocity of particle number change in 4m-dimensional two-particle elementary 
phase volume. This change is caused by the interactions between particles. Equation (1371) 
is the fractional generalization of the first Bogoliubov equation. If a = 1, then we have the 
first Bogoliubov equation for non-Hamiltonian systems [39] . 



4 Vlasov equation with fractional derivatives 

Let us consider the particles as statistical independent systems. Then 

P2(qi,Pi,q2,P2,t) = pi(qi,pi,t)pi(q2,P2,t)- (39) 
Substitution of (l39l) into ( l38l) gives 

/(P2) = -D;^^pi/[2]Fi2Pi(q2,P2,t), (40) 



where pi = pi(qi,pi,t). 
Let us define 

Then, 



F^^^(qi,Pi,t) =/[2]Fi2Pi(q2,P2,t). 



/(P2) = -D^,(piF^^0. (41) 
Substituting of (HT]) into (137|) . we obtain 

^ + D^^(VipO + D^^ ((F^ + {N- 1)F^/0pi) = (42) 

that is a closed equation for the one-particle distribution function with the external force 
F^ and the effective force F^-^'^. Equation 0421) is the fractional generalization of the Vlasov 
equation [3ll[35] that has phase-space derivatives of noninteger order. For a = 1, we get the 
Vlasov equation for the non-Hamiltonian systems |39j . 
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5 Fokker-Planck-Zaslavsky equation for phase-space 

The Fokker-Planck equations with fractional coordinate derivatives have been suggested by 
Zaslavsky [12] to describe chaotic dynamics. It is known that Fokker-Planck equation can 
be derived from the Liouville equation |3nillllll2]- In this section, we obtain Fokker-Planck- 
Zaslavsky equation that has fractional derivatives in phase space. 

Let us consider a system of N identical particles and the Brownian particle that is 
described by the distribution function 

pn+1 = P7v+i(q,p, Q,^',^), 

where 

p = (pi, ...,PAr), Pk = iPkl,-;Pkm) 

are the coordinates and momenta of the particles; Q = (Qs) and P = (Pg) {s = 1, ...,m) are 
Brownian particle coordinates and momenta. The normalization condition is 

i[l,...,N,N + l]pN+i = l. (43) 

The distribution function for the Brownian particle is defined by 

Pb{Q, P, t) = /[I, iV]p;v+i(q, P, Q, P, t). (44) 

The Liouville equation for Pn+i is 

^P^-^{LM + LB)pN+l = 0, (45) 

where 

N,m 

- ^LnP = Yl {^L (G'sP) + (^.V)) , (46) 



k,s 
N,m 

11 

k,s 

Here, and Lb are Liouville operators with fractional derivatives, and 



iLbP = J2{^t(9sP) + ^%{fsP))- (47) 



B'XB = (D'XA)-^ D'XB. 
The functions and are defined by the equations of motion for particle, 

^ = Gj(q,p), ^=F^{ci,p,Q,P), k = l,...,N. (4J 

9 



The Hamilton equations for the Brownian particle 



^ = gs{Q,P). ^ = fs{ci,p,Q,P). (49) 

define gs and fs- 

Let us use the boundary condition in the form 

lim pAr+i(q,p,Q,P,t) = pAr(q,p,Q,T)pij(Q,P,t), (50) 

where 

p^(q, p, Q, T) = exp f3{T - Hiq, p, Q)) (51) 
is the canonical Gibbs distribution for 

AT 

H{q, p, Q) = H^iq, p) + J] f/B(qfc, Q)- (52) 

fc=i 

Here, Hn is a Hamiltonian of an n-particle system, and Ub is an energy of interaction 
between particles and Brownian particle. If we suppose 

= pks/m, gs = Ps/M, (53) 

then 

N,m 2 

Hn{ci, P) = J2i^ + Yl ^(q'^' ^')- (54) 

k,s k<l 

The boundary condition (l50l) can be realized [36j by the infinitesimal source term in the 
Liouville equation: 

-^-^^ - i{LN + Lb)pn+i = -e{pN+i - PnPb)- (55) 
Integrating (!55|) by I[1,...,N], we obtain 

^ + J] D^^(^,ps) + J[l, AT] J] D?;^(/,p^+i) = 0, (56) 

which is the Liouville equation for reduced distribution function of the Brownian particle. 
The formal solution of Eq. fl55|) has the form 



pN+l{t) = PB{t)pN 



f dr e'r^-^r(L^+Ls) _ + LB))pB{t + t)pn. (57) 
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Substituting (J57l) into (l56i) . we get 

r^ rrt m 

s=l s=l 

-r[l,...,N]Y,^% / rfre^^e-^^(^-+^^)(— -z(LA. + LB))pB(t + r)pjv = 0. (58) 

1 J —oo \OT 



s=l 



The expression J[l, .., A^]/sPAr can be considered as average value of the force fs- For the 
canonical Gibbs distribution flSTl) it is equal to zero. Using 

D^^p^v = ^fi'^PN: (59) 
where /i^'' is a fractional potential force [25 • 

= -D^^f/B, (60) 

we have 



HbPn+i = ( -JTTTfPB + Dg^ (^,pb) + B'^^ifsps) ) PiV- 



It can be proved by integration that the term 

dpB 



dt 

in Eq. fl58l) does not contribute. Then fl58l) gives 



+ Dq.(^?.Pb) (61) 



^P£ 

at 



//t //t /.(J 

+ 5^D^^((7,PB) + 5^D^/[l,...,iV] / rfre-/,e-' 



iriLN+Ls) p^. 



■ (D^^,(/.,pB(t + r)) + ^^PB(t + r)) =0. (62) 

Equation ( |62l) is a closed integro-differential equation for the distribution function ps- Note 
that fs can be presented as 

where /i^"* is a potential force (1^0]) . and /i"^ is a non-potential force that acts on the Brownian 
particle. For the equilibrium approximation P ~ {MkTy^'^, iLs ^ M~^/^ and iL^ ^ m~^/^. 
If M >> m, we can use perturbation theory. 

Using the approximation psit + t) = psit) for Eq. ( l62ll . we obtain 

^ + E ^Qs(9sPb) + Yl Dp. (^D^^,(7L,pB(t)) + llPs'PBit)) = 0, (63) 

s=l s=l 
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where ^ 

7],, = /3Mr [1, N] I dr e'-fse-'^'^^fs'PN, (64) 

J —00 



7L, = PMP[1, ...,N] f dr e-/,e-^-^-/irV^- (65) 



00 



If /, = /i^), then 

lis' = lis'- As a resuh, we derive the Fokker- Planck- Zaslavsky equation 
, [H] for the phase space. 



6 Linear fractional kinetic equation 

Let us consider Eq. (1371) with /(P2) = 0, Vi = p/m = v, and F*^ = eE, B = 0. Then 

^ + (v,D^pO+e(E,D>i) = 0, (66) 

where 

m 

(v,D>)=5^(i;„D^^p). (67) 

s=l 

If we take into account the magnetic field (B 7^ 0), then we must use the fractional 
generalization of Leibnitz rules 

Dg(/^) = E r(.. + i)r(a-.. + i) <P'^"^)°'^^- 

s=0 \ / \ / 

where s are integer numbers. In this case, Eq. fl66l) has the addition term 

([p, B]pO = ^ D° {e,,^p,B^,p,) = ^J2 ^^iraB^^K^ {p,p,) = 

klm Mm 

mc^-^ i u + 1 i a — ? + 1 

klm i=0 \ / \ / 



e 

mc 



^£kimB,n {[^p^pi]pi + a[D^^ V]40 = 

fc/m 

— $^e.^Ap^[D^,Pi] = — ((D^^pi), [p,B]) . (69) 

klm 

Let us consider the perturbation [23 [M] of the distribution function in the form 

pi = pi + 6pi{t,q,p), (70) 
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where pi is a homogeneous stationary density of probabihty that satisfies Eq. (]66l) for E = 0. 
Substituting ([70]) into Eq. ([66]), we get 

^ + (v, D^Vi) + e(E, D>~i) = 0. (71) 

Equation (1711) is hnear fractional kinetic equation for the first perturbation 6pi of the distri- 
bution function. Solutions of fractional linear kinetic equations ( 17T1) are considered in Ref. 
[13]. For E = 0, the function 6pi is described by the function 



(^7,t)-V"L„ [q,{g,t)-'/^] , (72) 

where = Vs{D'^^qs)'^ , and 

1 r+oo 

L^[x\ = — dk e-^'="e-"l'=l" (73) 

is the Levy stable p.d.f. |13]. The examples of Lj^x] are shown in Fig. 1. 
For a = 1, the function ([73]) gives the Cauchy distribution 

Li[x] = -^^, (74) 

TT + 1 



and ([72]) is 

1 {9st)-^ 



For a = 2, Eq. ([75]) gives the Gauss distribution: 

1 



and the function ([72]) is 



Tlx — ' n\ 

n=l 



(75) 



UV\ = (76) 
(^,t)-V2j-e-'^^^/(^^^*). (77) 

2y/TT 

For 1 < a < 2, the function La[x] can be presented as the expansion 

The asymptotic (x — oo, 1 < a < 2) is given by 

LJx] y^(_i)«a;-"£il±!^ sin(n7r/2). (79) 

nx ^-^ n\ 

n=\ 



- > (-x)" ^ / ^ sin(n7r/2). (75 



As the result, the asymptotic of the solution, exhibits the power-like tails for a; oo. This 
tails is the important property of solutions of equations with the noninteger derivative. 
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Figure 1: Gauss p.d.f. {a — 2), Levy p.d.f. {a — 1.6), and Cauchy p.d.f. {a — 1.0). Levy 
for a — 1.6 lies between Cauchy and Gauss p.d.f. In the asymptotic x — > oo and x > 3 on 
the plot, the upper curve is Cauchy p.d.f., and the lower curve is the Gauss p.d.f. 
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7 Conclusion 



In this paper, we consider equations with derivatives of noninteger order that can be used 
in statistical mechanics and physical kinetics. We derive the Liouville, Bogoliubov, Vlasov 
and Fokker-Planck equations with fractional derivatives with respect to coordinates and 
momenta. To derive the fractional Liouville equation (FLE), we consider the conservation 
of probability in the fractional differential volume element. This element can be considered 
as a small part of the phase-space set with noninteger dimension. Using the FLE, we obtain 
a fractional generalization of the Bogoliubov hierarchy equations. These equations describe 
the evolution of the reduced density of probability in the fractional phase-space volume 
element. Fractional Bogoliubov equations can be used to derive fractional kinetic equations 
[5l [121 [13] • In this paper, we obtain Fokker-Planck-Zaslavsky equation, fractional Vlasov 
and linear kinetic equations. 

The fractional kinetics is related to the equations that contains derivatives of noninteger 
order. These equations appear in the description of chaotic dynamics, and the fractal media. 
The fractional derivatives can be connected with long-range power-law interaction of the 
systems |23l [MJ [25]. For noninteger derivatives with respect to coordinates, we have the 
power-like tails as the important property of the solutions of the fractional equations. 
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